Quantum time of flight distribution for cold trapped atoms 
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The time of flight distribution for a cloud of cold atoms falling freely under gravity is considered. 
We generalise the probability current density approach to calculate the quantum arrival time distri- 
bution for the mixed state describing the Maxwell-Boltzmann distribution of velocities for the falling 
atoms. We find an empirically testable difference between the time of flight distribution calculated 
using the quantum probability current and that obtained from a purely classical treatment which is 
usually employed in analysing time of flight measurements. The classical time of flight distribution 
matches with the quantum distribution in the large mass and high temperature limits. 
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I. INTRODUCTION 

In recent times laser cooling and trapping of atoms has 
become an area of active research [TJ. The measurement 
of the initial temperature of the cloud of atoms is crucial 
for characterising the properties of atom traps. The tem- 
perature of the cloud can be inferred from the velocity 
distribution of atoms in the cloud. A well-known tech- 
nique of measuring this velocity distribution is the time 
of flight (TOF) method. Measurements of the TOF dis- 
tribution have been employed to analyse various experi- 
mental data such as those involving ions and isotopes 0, 
and also inperforming mass spectroscopy of biomolccules 
like DNAgf. 

The theoretical treatment of the TOF distribution that 
can be obtained using for instance, the Green's func- 
tion method however, turns out to produce perfect 
agreement with the TOF distribution obtained by us- 
ing Newton's equations for ballistic motion of particles 
accelerated by the earth's gravitational field[j|. Thus, 
the interpretation of the results of the various TOF 
experiments (T), [H, Q where classical trajectories are in- 
ferred from Newtonian mechanics^] remains debatable, 
especially in the domain of small atomic masses and low 
temperatures where quantum mechanical effects should 
be significant. 

Though there exists no unique prescription for the def- 
inition of time of flight and arrival time in quantum me- 
chanics, experimentalists measure arrival times of ele- 
mentary particles, atoms and molecules using the TOF 
methods. In spite of the difficulties to give time an ob- 
servable status in quantum mechanics, several logically 
consistent schemes for the treatment of the arrival time 
distribution have been formulated, such as those based on 
axiomatic appr aches Q, opearator const ruct ions [H, and 



trajectory models[9)]. It is thus desirable that some of 
the conceptually sound theoretical formulations of the 
quantum mechanical arrival time distribution [Io| be con- 
fronted with accurate experimental data. If such quan- 
tum mechanical approaches are employed for analysing 
experiments using TOF measurements, it should not only 
enable to determine the empirical viability of various 
competing arrival time models (Toj|. but also possibly shed 
new light on the conventional interpretation of the results 
of these experiments. 

In this paper we employ the probability current 
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approach llj towards obtaining the quantum arrival time 
distribution of cold trapped atoms. The probability cur- 
rent approach for computation of the mean arrival time 
of a quantum ensemble not only provides an unambigu- 
ous definition of arrival time at the quantum mechanical 
levelpj], [ID, [l3j], but also adresses the issue of obtaining 
the proper classical limit of the time of flight of massive 
quantum particles [3, EH ■ Here we derive the quantum 
arrival time distribution for the case of initially trapped 
atomic clouds that are subsequently allowed to fall freely 
under gravity 0, EH- We compute the mean time of ar- 
rival for these atoms and compare it with the mean ar- 
rival time obtained through the classical time of flight 
analysis Q that has frequently been employed for such 
experiments [l|, 0, 3- Our analysis predicts the mass and 
temperature range of the atomic clouds where the quan- 
tum mechanical treatment alters the arrival time distri- 
bution and the mean arrival time from that obtained 
through the classical analysis. 



II. CLASSICAL ANALYSIS OF TIME OF 
FLIGHT MEASUREMENTS 

We begin with a brief description of the classical anal- 
ysis of TOF measurements of trapped atoms. A probe 
laser, focussed in the form of a sheet, is placed under- 
neath the atomic cloud. When the trapping forces are 
turned off, the cold atom cloud falls through the laser 
probe under the influence of gravity. It is then possi- 
ble to detect the fluorescence from the atoms as they 
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reach the sheet. The fluorescence is measured as a func- 
tion of time and the initial temperature of the cloud is 
determined by fitting the experimental result to the the- 
oretically predicted TOF signal of the cloud [J, |T(|. A 
detailed derivation of the TOF signal recorded by the 
detector (that is, the number of atoms arriving at the 
probe laser as a function of time) was derived by Yavin 
et a/0. 

The cloud of atoms consisting of noninteracting parti- 
cles has a Maxwell-Boltzmann velocity distribution given 
(in one dimension) by 



U(v)dv = ( 
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where T is the initial temperature of the cloud, and m 
is the atomic mass. Using the Newton's equations for 
ballistic motion of a particle accelerated by the earth's 
gravitational field (in the vertical z-dircction), the veloc- 
ity is obtained in terms of the time of flight as 



v = (z + l -gt 2 )/t 



(2) 



Substituting the above expression for v from Eq.(2) in 
Eq.(l), one can obtain the time of flight distribution at 
an arbitrary distance z, given by 
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The corresponding classical mean time of flight or mean 
arrival time tc for the atomic cloud calculated using 
lie (£) as the time of flight distribution is given by 



f~Ilc(t)tdt 



(4) 



For simplicity, here we restrict ourselves to the case of a 
point-sized cloud. The three-dimensional calculation was 
done by Yavin et al[f| using a simple coordinate trans- 
formation and the same expression was obtained for the 
TOF distribution. These authors Q have claimed perfect 
agreement of their results with a previous calculation [J| 
where the TOF distribution was derived using a sophisti- 
cated Green's function technique. We call this TOF dis- 
tribution given by Eq. ^ as the classical time of flight dis- 
tribution, and Eq. (fJJ denotes the corresponding classical 
mean arrival time since the classical Newtonian equation 
is used from the outset to derive this TOF distribution. 



III. THE QUANTUM ARRIVAL TIME 
DISTRIBUTION THROUGH THE 
PROBABILITY CURRENT 

Our aim here is to derive an expression for the time 
of flight distribution for the atomic cloud through the 



quantum probability current without using any classical 
ingredients. To that end, let the initial state of each of 
the atoms be represented by a one dimensional Gaussian 
wave function of the form 
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centered at z = and moving with a group velocity v. 
The Schrddinger time evolved wave function under the 
Hamiltonian H = p 2 /2m + mgz is given by 
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where St — (Jq (l + ihtj2raa 2 ^. 

Considering the free fall of the atoms under gravity, the 
expression for the Schrddinger probability current density 
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for the time evolved state is calculated using the initial 
state given by Eq.(f5]) to be 



J{z,t) = P(z,t) 
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where the expression for the position probability distri- 
bution is given by 
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with a = a (l + h 2 t 2 /4m 2 a%) 1/2 . The modulus of the 
probability current density J{z,t) given by Eq.® pro- 
vides the arrival time distribution for a pure wave packet 
falling under gravity. Note that the quantum probabil- 
ity current as defined by Eq.([7]) is formally ambiguous 
up to a total divergence termpj}. However, J(z,t) can 
be uniquely defined through relativistic wave equations 
which impart appropriate spin-dependent corrections to 
it that persist even in the non-relativistic limit (l2l. [l8| . 
The ensuing arrival time distribution defined through the 
probability current thus contains a spin-dependent cor- 
rection for particles with spin [3]. 

The atomic cloud is represented by an ensemble of par- 
ticles in thermal equilibrium with a thermal distribution 
of initial velocities. Each particle has a wave function of 
the form ([5]), with a Maxwell-Boltzmann distribution of 
initial velocities given by Eq.(l). Thus the initial thermal 
state of the atomic cloud we have described is a mixed 
state. We obtain the corresponding position probabil- 
ity distribution by averaging the pure state distribution 
© over a thermal distribution of initial velocities. The 
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FIG. 1: The classical (U c (t)) and quantum (H Q (t)) TOF dis- 
tributions of the atomic cloud falling freely under gravity are 
plotted for varying mass of the atoms at a fixed temperature 
T — 2.5 x 10~ 6 K with a = l(T 5 cm, Z = -30cm. 



FIG. 2: The classical (Ilc(i)) and quantum (n Q {t)) TOF 
distributions of the atomic cloud falling freely under gravity 
are plotted for varying temperatures at a fixed mass of Rb 
atom (m = 85. 4678a. m.u) with ao = lCP 5 cm, Z = —30cm. 



mal distribution of initial velocities. The result is 
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where or 2 = cr 2 + (kT/m)t 2 . The peak of the position 
probability distribution Pt(z, t) follows the classical tra- 
jectory and the effect of the mass and temperature de- 
pendences of the position probability occurs essentially 
because the spreading of the wave packet is different for 
different atomic mass and temperature of the cloud. 

Taking the modulus of the quantum probability cur- 
rent density as determining the quantum arrival time 
distribution [ 11], we obtain the arrival time probability 
distribution for the atomic cloud (which is a mixed state) 
by averaging the pure state distribution ((5|) over a ther- 
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In this way we generalise the probability current density 
approach to calculate the quantum arrival time distribu- 
tion for the mixed state at finite temperature. It may be 
mentioned here that in the present calculation we neglect 
the small spin-dependent correction [l3j that may appear 
in the probability current density, as mentioned above, 
and consequently the mean arrival time that we compute 
for any fermionic atoms. 

The corresponding quantum mean arrival time tq 
which is an observable quantity, is given by 



_ = / °°n Q (t)trft 
TQ J Q °°n Q (t)dt 



(12) 



One may note that though the integral in the numera- 
tor of Eq. (fl"2")) formally diverges, several techniques have 
been employed in the literature ensuring rapid fall off 
for the probability distributions asymptotically so 
that convergent results are obtained for the integrated 
arrival time. For our present purposes it is sufficient to 
employ the simple strategy of taking a cut-off (t = tc) 
in the upper limit of the time integral [3, O with 
tc = y/2(Z + 6<7f C ) / g where ate — ^r{t = tc) is the width 
of the wave packet at time tc. Thus, our computations 
of the mean arrival times are valid up to the 6cr level of 
spread in the wave function. 

If we impose now the classical limit for the quan- 
tum TOF distribution given by Eq. lfTTj) . then one can 
check that under the large mass and high tempera- 
ture limits, and when do << (kT/m)t 2 , one can take 
a T 2 « {kT/m)t 2 . Thus, U Q (t) = U c (t), i.e., the 
two distributions match in the limit of large mass and 
high temperature. The probability current method of 
computing the quantum arrival time distribution fur- 
nishes an effective way of approaching the classical limit 
of the distribution by smoothly varying the parame- 
ters such as mass and temperature of the quantum 
distribution [3, EH]. Note also, that the mass depen- 
dence of arrival time distribution given by Eg. (|1 1[) and 
consequently, the observable mean arrival time given by 
Eq. (|12|) . signifies the quantum mechanical violation of 
the gravitational weak equivalence principle (l5j. Thus 
TOF measurements [1 0, [|, H [lj] offer a practical possi- 
bility for experimental demonstration of the equivalence 
principle violation at the quantum level [2fj|. 



IV. NUMERICAL RESULTS 

We perform a numerical study of the quantum arrival 
time distribution of the falling atomic cloud by varying 
its mass and temperature separately. We first plot lie (t) 
and ng(i) for fixed temperature of 2.5 x 10~ 6 K in Fig.l. 
It is seen the classical and the quantum distributions are 
clearly different for clouds of small atomic mass such as 
Li and Na. However, as one increases the atomic mass, 
one sees that lie (t) and Uq (t) begin to overlap for heavy 
atoms such as Rb for this value of temperature. The 
temperature variation of the arrival time distributions are 
displayed in Fig. 2 where one sees that even for heavy Rb 
atoms, the two distributions are quite distinct in the low 
(nano kelvin) temperature range. It would be interesting 
if our prediction in the low temperature and lower atomic 
mass region where the quantum time of flight distribution 
sharply differs from the classical TOF distribution could 
be verified in actual experiments. 

The variation with mass of the quantum and classical 
mean arrival times at a particular detector location for 
the ensemble of falling atoms is depicted in Fig. 3. One 
can see that in the limit of large mass the mean arrival 
time tq asymptotically approaches the classical result. 
One can also investigate the variation of the mean ar- 
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FIG. 3: The mass variation of the mean arrival times tq and 
tc are shown in the figure for a fixed value of temperature 
(T = 1.41 x 10" 6 iO with Z = -30 cm. The quantum mean 
arrival time tq is plotted for two different values of <tq. 2(a) 
tq for a = lCT 5 em, 2(b) tq for a = 2 x 10" 6 cm, 2(c) The 
mean arrival time tc calculated through the classical TOF 
distribution Uc(t) with T = 1.41 x lO" 6 ^ and Z = -30 cm. 

rival times (tq and Tc) with varying temperature of the 
cloud, and obtain similar results, as expected from the 
temperature variation of the classical and quantum ar- 
rival time distributions plotted in Fig. 2. The mass and 
temperature dependences of the quantum arrival time 
distribution and consequently the quantum mean arrival 
time arise essentially due to the spread of the wave packet 
for the atoms. A smaller value of the initial width <ro for 
the wave packet results in its faster spread. The amount 
of departure of the quantum distribution from its clas- 
sical counterpart is thus contingent on the magnitude 
of the ensemble spread (since Hc(t) is independent of 
(Tq). This is clearly depicted in Fig. 4 where the classi- 
cal (Uc(t)) and quantum (ilg(t)) TOF distributions are 
plotted for three different values of <tq at a fixed mass 
and temperature. 



V. SUMMARY AND CONCLUSIONS 

To summarize, in this work we have considered the 
analysis of the time of flight measurements of falling cold 
atomic clouds. The inference of the temperature of the 
cloud in various experiments [l], 0, El is usually performed 
through a classical analysis in which the results obtained 
are same as through the solution of Newton's equations 
for ballistic motion of particles falling under gravity Q. 
Here we emphasize the relevance of employing a quan- 
tum mechanical arrival time distribution for the analysis 
of such experiments. We use the probability current den- 



5 



(a) 




■6 

(J. - 10 cm 






->IIc 











(b) 




i 

O, = 5xl0 « 

- 


L 


J 






A i < i 



a 

| 20- 



(0 


->Hc 




W 




o, ■ io' 5 » 1 






n 


Cj = S i' io" 5 ™ 


n Q >H 






i . 










6 

| 20 








L 


u 




' 


) 



O.i 0.15 2 25 O.J 0.35 0.4 



0.15 0.2 0.25 0.3 0.35 0.4 



Time :ri sec 



lire live;. 



FIG. 4: The classical (U c (t)) and quantum (Jl Q (t)) TOF 
distributions of the atomic cloud falling freely under gravity 
are plotted for four different values of <tq at a fixed mass of 
Be atom (m = Q.Ola.m.u) and at a fixed temperature T = 
3.0 x 10" 6 ^ with Z = -30cm. 



sity approach towards obtaining the arrival time or the 
TOF distribution. Our definition of the quantum ar- 
rival time distribution and the observable mean arrival 
time in terms of the modulus of the probability current 
density is particularly motivated from the equation of 
continuity, and other physical considerations discussed 



in the literature d [TT|, [H, d, H- Further, we gener- 
alise the probability current density approach to calcu- 
late the quantum arrival time distribution for a mixed 
state describing the Maxwell-Boltzmann distribution of 
velocities for the falling atomic clouds in the relevant 
experiments [ll, 0, EH- We compute the TOF distribu- 
tion and mean arrival time through this scheme and com- 
pare our results with those obtained through a classical 
analysis @ for various atomic masses. 

The obtained quantum arrival time distribution 
matches with the classical TOF distribution in the high 
temperature and the large mass limits, hence furnishing 
another example of the smooth emergence of the classical 
limit of the quantum arrival time in the framework 
of the probability current approach. However, a clear 
distinction between the quantum and the classical distri- 
butions is exhibited for either small atomic mass, or low 
temperature of the cloud. This results from differential 
wave packet spreading depending upon the mass, velocity 
and width of the wave packet for the atoms. Our scheme 
thus provides a new method for experimental verification 
of the probability current density approach for calculat- 
ing the arrival time distribution, in addition to an earlier 
proposed method using the spin rotator as a quantum 
clock[22|. Finally, we wish to emphasize that more inves- 
tigations of modern experiments employing time of flight 
techniques should be performed using various quantum 
mechanical schemes [10J. Such studies have the potential 
to empirically resolve ambiguities inherent in the theoret- 
ical formulations of the quantum arrival time distribution 
using cold trapped atom experimental techniques, and 
may also shed new light on the inference of the experi- 
mental data. 
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